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SUMMARY 


Problems relating to the influence of lunar and solar per- 
turbations on the motion of artificial satellites are analyzed by 
an extension of Cayley’s development of the perturbative function 
in the lunar theory. In addition, the results are modified for 
incorporation into the Hansen-type theory used by the NASA Space 
Computing Center. The theory is applied to the orbits of the 
Vanguard I and Explorer VI satellites, and the results of detailed 
computations for these satellites are given together with a physi- 
cal description of the perturbations in terms of resonance effects. 
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DEVELOPMENT OF THE LUNAR AND SOLAR 
PERTURBATIONS IN THE MOTION 
OF AN ARTIFICIAL SATELLITE* 


INTRODUCTION 


This paper is concerned primarily with problems related to the influence of lunar 
and solar perturbations on the motion of artificial satellites. The importance of these 
problems is indicated by Kozai's discovery that the perigee height and lifetime of a sat- 
ellite may be strongly affected by these perturbations. 

The basis for the computations is provided by an analytical development of the dis- 
turbing function, which is an extension of Cayley’s development of the solar perturbative 
function in the lunar theory. The relations between perigee -height variations and launch 
conditions have been investigated by using a modification of this development. Values of 
perturbations in the perigee height for the satellites Vanguard I (1958 fi^) and Explorer VI 
(1959 h) were computed from the resultant trigonometric series. A program was devel- 
oped for computing lunar and solar perturbations with the aid of the IBM 704. This pro- 
gram permits the inclusion of any value of the eccentricity and of the inclination, and gives 
the variation of the perigee height in the form of a trigonometric series with numerical 
coefficients. Also, the perturbations in Hansen's coordinates have been computed; this 
program can be included in the existing Vanguard scheme for general oblateness per- 
turbations. The formulas used in these computations can be applied to the development 
of a Hansen-type numerical theory for an artificial satellite. 


DEVELOPMENT OF THE DISTURBING FUNCTION 

The importance of determining the lunar and solar perturbations in the motion of an 
artificial satellite was indicated by Kozai's discovery that certain long-period terms in 


^Originally presented at the First International Space Symposium, sponsored by the Com 
mittee on Space Research (COSPAR), Nice, France. January I960. 
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the development of the disturbing function cause large perturbations in the elements and 
that, in this way, the orbital lifetime of the satellite can be considerably affected (Refer- 
ence 1). In the present treatment, the analytical development of the two main terms of 
the disturbing function are expressed 






)■ 


where 


s = cos ( r, r ' ) . 


In this expression ( r, r') is the angle between the position vectors of the satellite and 
of the perturbing body, either the sun or the moon, with respect to the earth. The main 
term is the second Legendre polynomial in the harmonic expansion of the disturbing 
function. The third Legendre polynomial is known as the parallactic term. 

The analytical development is important because it permits investigation of the per- 
turbations for an entire group of satellites having similar elements and, in particular, is 
useful in investigating problems connected with the effects of resonance, i.e., stability 
problems. 


The problem of developing the disturbing effects caused by the sun and the moon is 
related to the lunar problem but with the difference that the orbital inclinations of both 
the disturbed and the disturbing body to the basic reference plane can be large in the 
case considered here. The arrangement given for the development of the disturbing 
function, based on Cayley's work (Reference 2), permits the inclusion of any power of the 
eccentricity and is valid for all inclinations (Reference 3). 


For the time being, only long-period terms are evaluated, and the terms depending 
on the mean anomaly of the satellite are excluded. The resultant "abbreviated version" 
is given on pages 25 through 27. This last development was used to investigate the influ- 
ence of the combined effect of drag and lunar and solar perturbations on the orbital life- 
time of the satellite and in further investigations of the orbits of satellites in the NASA 
space research program. Future satellite orbits undoubtedly will have elements such 
that the terms depending on the mean anomaly of the satellite and the higher powers of 
the eccentricity can be expected to become more important, especially regarding the 

possible development of resonances associated with the commensurability of satellite 
and lunar periods. 
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The perturbations may also be developed by a purely numerical method based on the 
use of fast computing machines. In this method the values of the orbit elements a , e , 
i, a' , e' , i' are substituted into the coefficients. The numerical method is convenient 
in developing the perturbations according to Hansen ^s theory and in investigating the 
variation of the perigee distance. A program for the development in terms of the eccen- 
tric anomaly and a program for the development in terms of the mean anomaly have been 
prepared. Both programs can help to supply information concerning the optimum condi- 
tion for launching in connection with the solar and lunar perturbations. 

Direct numerical integration of the equations of motion with the solar and lunar per- 
turbations, by means of a computer program, has confirmed the results of the anal5dical 
treatment and the Fourier series development. 


COMPUTATION OF SATELLITE PERTURBATIONS 


The following notations are used: 


a 

e 

i 

Q 

g 

r 

X, y, z 
r 
f 

r 


semimajor axis of satellite orbit 
eccentricity 

angle of inclination to equatorial plane 

argument of perigee 

right ascension of ascending node 

mean anomaly 

position vector 

rectangular coordinates 

radius 

true anomaly 
i 

sin . 


The corresponding elements of the disturbing body are designated by primes; a ' , e ' , 
etc. The equatorial plane is taken as the basic reference plane, and the precession and 
nutation of the earth's axis are neglected. The disturbing function has the form 


n 




^ P^CS) + ^ P3(S) 


( 1 ) 


or 


n 


( 2 ) 
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In the latter expression 


S = cos 


(r, r') (t + co + o - f ' - oj' ) 


+ 7^ (1 “7^^) cos (f + a) - 0 + +6j') 


+ y'^ (1 “ 7^) cos (f+cj + ^ + f +cu) 


+ 7 ^ 7 ^ ^ 


cos ( f + 


f ' - cj') 


2yy' yjl - 7^ x/l - 7^2 cos (f + oj - f - ) 


2 77' \/l " 7 ^ 


— 7^2 cos (f +^x>+f^ +6c»*) 


where 


6> = n - n' 


and P2, P3, 
considered: 


are Legendre polynomials. Only the first two terms of Equation 2 are 




If i ' / 0, then it is convenient to represent Equation 4 in the form of Radau (Reference 4 ): 

„ m'a^/rV/a^\^/l r^- 3 z^ 


\s) \r^) \4 


3 X ^ ~ y ^ xy x y xz x z 

4 2 I 2 ^ 2 f 2 ^ 2 > 2 


yz \ 
r2 r'2y 


m'a^/r\^/a '\^/3 r^x“5xz^ r'^x'“5x'z'^ 

TTT U j j \ "8 “3 


3 r^y - Syz^ r'^^y' - Sy'z'^ 1 3 r^z - Sz*^ 3 r'^z' - 5 z"' 
" 8 + 4 7 ’ 773 


^ - 3 xy 2 x'^ - 3 x'y'' 

* 8 73 


t s 3x2y - y3 3x'2y' - y'3 
8 ,3 ,.3 


15 x2z “ v^z x'2z' - V^^z 


4 ^3 


xyz X ' y ' z ' \ 

7 ^ ) ’ 


( 5 ) 
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and to substitute 

— zz (1 - cos (f + O) + ^) + 7^ cos (f + w - , 

Z = (1 - 7 ^) sin (f + 60 + ( 9 ) - -y^ sin (f + , 

Z = 2 7 y/l - 7 ^ sin ( f + a;) , 


Z_ = cos ( f ' + ^* ) ) 
r 

/ 

= (1 “ 27 '^) sin + 60' ) , 

_ 5 _ = 27' y 1 - 7' 2 sin (f + ) • 


Substitution of these expressions into Equation 5 results in the following development for 
n in terms of the true anomalies: 


afi 



(sum of terms of Types I-V) 


(sum of terms of Types VI-XIII): 


( 6 ) 


Type I 

+ -i (1 - 67 ^ + er**) (1 - 67 '^ + 6r"*) 

+ 3 ( 7 ^ - 7^) ( 7 '^ ” 7 '^) cos 26? 

+ 3 (7 “ 27 ^) (7* ' 27 '^) \/l - ^ \/l - 7 '^ cos 0 ; 

Type II 

+ |_ (y2 _ ^ 4 ^ ^^'2 _ cos (2f + 26*J - 2f ' - 2co') 

+ A (1 - y2^2 _ 7'2)2 ^2f + 2 (^ “ 2f ' - 2661^ + 20 ) 

+ -^7^ 7 "^ cos (2f + 260 “ 2f ^ “ 26 o' - 2/9) 

+ 3 (7 - 7^) ( 7 ' " 7 '^) yj\ ~ 7 '^ y 1 - 7 ' ^ cos (2f + 2^ ” 2f - 2co' + 


+ 37 ^ 7 '^ y 1 - 7 ^ y 1 - 7^2 cos (2f + 264^ ^ 2f^ "• 2co' - 9 ) ; 
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Type lU 

+ |- (r^ - 7 *) (1 - 67' 2 + 67"*) cos (2f + 2 co ) 

+ y (1 - 7^)^ (7'^ - y"*) cos (2f + 2w + 20) 

+ |- 7 ^ (7'^ - 7 "*) cos (2f + 2a. - 20) 

- 3 (7 - 7^) (7' - 27'^) \/l - 7^ \/l - 7' 2 cos (2f + 2o) f 0 ) 

+ 37^(7' - 27'^) v/l - 7 ^ yi - 7'^ cos ( 2 f + 2 o) - 0 ) ; 

Type IV 

+ y (7'^ - r"*) (1 - 67^ + 67'*) cos (2f' + 2a.') 

+ ^( 7 ^ - 7“) (1 - 7'^)^ cos (2f ' + 2a.' - 20) 

+ y (7^ - 7“*) 7'^ cos (2f + 2 o)' + 20 ) 

+ 3 (7 - 27^) 7'^ \/l - 7^ \/l - 7'^ cos (2f' + 2a.' + 0 ) 

- 3 (7 - 27^) (7' - 7'^) ■JX - 7^ 7 l - 7'2 cos (2f + 2a' - 0) ; 

ry/>c y 

+ y( 7 ^ - 7'*) (7'^ - 7"*) cos (2f + 2a. + 2f' + 2a.') 

3 2 

+ “ 7^) 7"* cos (2f + 2oj + 2f' + 2oj' + 20) 

+ -f- 7 '* (1 - 7'^)^ cos (2f + 2o) + 2f' + 2o)' - 20) 

- 3 (7 - 7^) 7' 3 \}x - y'^ X - 7 ' 2 cos (2f + 2a. + 2f' + 2a>' + 0) 

- 37^(7' - 7’^) x/l - 7 ^ Vl - 7'^ 


cos (2f + 2oj + 2f' + 2oj' - 0) : 
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Type VI 

“*y(X “ 5/^ + 57^) {y* “ 5y'^ + 57'^)\/l - 7^ yj\ - 7' ^ cos (f +oj+f' +cj') 

+ |-(1 - 7^) (1 - IO72 + 157^) (67^2 - 207'"^ + 157'^) cos (f + o) + f* + + < 9 ) 

o 

+ -|-(6y^ ~ 20/'* + 15 y®) (1 - y'^) (1 - lOy'^ + 15 y"*) cos (f + + f' + o>' - 0 ) 

O 

+ ^ (7 ” 47^ + 37^) (27'^ ~ 37' \/l “7^ y 1 _ y' 2 + 2^?) 

+ ^ (27^ “ 37^) (7' -* 47'^ + 37'^) \/l “ 7 ^\/l -7'^ cos (f + o) + f' +oj' - 2 ^^) 

+ ^ (7^ - 27"* 7^) (7^"^ - 7'^) cos (f + oj + f ' + oj' + 3 ^?) 

o 

+ ^ (r"* - r^) (r‘^ - 27"* + r'®) cos (f + r j + f' + <y - 30 ) ; 

o 

Type VII 

+ -^ ( y - 5y^ + 5y®) (y' - 5y' ^ + 5y'^) \/ 1 — y^ yj I - y' 2 (.Qg - -oj') 

+ J-(l - r^) (1 - 10y2 + 15y'‘) (1 - r'2) (1 - 10/'2 + 15y'^) cos ({ * a> - f - co' + 0) 

+ 1 ^( 67 ^ - 2o 7 + 15y®) (6y'2 - 20y"* + ISy'®) cos ( f + a; - f ' - oj' - 0) 

O 

+ ^ (t - 4y3 + 3yS) - 4 y'^ + 3 y'®) x/l - y^ cos (f + oj - f ' - w' + 20 ) 

+ ^ ( 2 y^ - 3 y®) ( 2 y'^ - 3 y'^) /l - y^ - y' 2 cos (f + - f' - w' - 20 ) 

+ ^ (r^ - 27*^ + y®) (y '2 - 2 y"* + y'®) cos (f + o; - f' - «' + 30 ) 

+ ^ (y'* - y^) (y"* - y'®) cos (f + cj - f' - <u' - 30 ).; 

Type VUl 

- ^ (y - 5y^ + 5y^) (y'^ - y'®) yj\ - y'^ \/l - y' ^ cos ( f + &j + 3f' + 3«') 

+ -| (1 - y 2 ) (1 - I 0 y 2 + ISy^) (y"^ - y'«) cos ( f + + 3 f' + 3 o)' + 0 ) 

+ ^ (6y^ - 20 y 4 + 15 y®) (y' - y'^)^ cos (f + + 3 f' + 3 &>' - 0 ) 

+ ^ (r - 4 y^ + 3 y^) y'^ -Jl - y'^ - y'^ cos (f + tu + 3f' + 3<i>' + 20) 

- ^ ( 2 y^ - 3 y^) y'(l - y'^)^ 7 1 - y^ ^/^ - y' ^ cos (f + o + 3 f' + 3w' - 20 ) 

15 2 

+ 3- (y - y^) y'® cos (f + oj + 3 f' + 3 ^' + 30 ) 

+ (y'* - y®) (1 - y'^)^ cos (f + w + 3f' + 3<>)' - 30 ) ; 



8 


Type IX 

+ ^ (r - 5 y 3 + SyS) (r'^ - 7 '®) v/l " 7 ^ V 1 " 7 ' ^ cos (f + a; - 3 f' - 3 a;') 

+ 11 (1 - -y 2 ) (1 _ 10^2 + 15 ^ 4 ) (^y> - y' 3)2 COS (f + a; - 3 f' " 3ai' + 0 ) 

8 

+ 11 (6y2 _ 207“* + 157®) (7'^ " 7'®) cos (f + a; - 3 f' - 3 a;' - B) 

8 

_ (7 - 473 + 375) (7' - 27'^ + 7'®) - y'^ Vl - 7'^ cos (f + - 3 f' - 3 ca' + 26 ) 

+ ^ ( 2 y^ - 3 r^) 7'^ yj\ - V 1 “ 7' ^ cos (f + c*> - 3 f' “ “ 2 ^) 

4 - 7^(1 “ 7^)^ (1 ” 7*^)^ + Cij ~ 3 f' ' 3 oi^ + 3 ^) 

8 

+ - 1 ^ (7^ “ 7^) 7 ^^ cos (f + a; ” 3 f * ” Sco *■ 3^-0 ; 

8 

Type X 

^ i|. (-^3 _ ^5) (■y* ^ 5^'3 + 5 ^^ 5^ Vl - 7 ^ yi - 7' 2 ( 3 f + 3 ^; + + a;' ) 

+ ^ 207'"^ + 157'^) cos ( 3 f + 3 c^ + f' + + 6 ^) 

8 

+ /^4 _ ^6^ _ 7'^) (1 ” 107'^ + 157'^) cos ( 3 f + 3 oj + f' + oi' - 0) 

8 

_ ^ (y - 27^ + 7 ^) (27'^ “ 37'^) yfl - 7^ \/l -7'^ cos ( 3 f + 3 oi + f' + O)' + 26 ^) 

+ ^ 7^(7' - 47'^ + 37'^) >/l “ 7 ^ \/ 1 -7'^ cos ( 3 f + 3 ^ 2 ) + f' + co' - 26 ^) 

+ (1 - 7^)^ (7'^ “ 7'^) cos ( 3 f + 3 ^^iJ + f' + + 36 *) 

8 

4- AA 7^^7^ ” 7* cos ( 3 f + 3 co + f ^ + C2J *" 3 ^) \ 

8 

Type XI 

+ (y3 _ ^5^ ^yf ^ ^yt3 4 . y/\ - 7^ yf^T- 7^ COS 3f + 3a; ” f' - o>' ) 

+ H (y - 7^)^ (1 “7'^) (1 “ 107"^ + 157'^) COS ( 3 f + 3 oj - f’ - +6^) 

8 

+ lA (y 4 _ y6^ (67'^ - 207'"^ + 15 /'^) cos ( 3 f + 3 a; - r - a;' - 6 *) 

8 

- ^ (y - 27^ + 7^) (7' "* 47'^ + 37'^) y 1 “ 7^ \/ 1 “ 7' ^ cos ( 3 f + 3 a) - f ' “ oj' + 2f>) 

+ ^ 7^(27'^ - 37'^) ~ 7 ^ Vl - 7' 2 cos ( 3 f + 3 o) - f' - o)' - 26 ?) 

+ (1 “ 7^)^ (7^ - 7'^)^ cos ( 3 f + 3^4; “ f' “ (2;' + 3 ^^) 

8 

+ — 7^(7^^ — 7*^) cos ( 3 f + 3o) “ f ' - (ii' - 36 ?) ; 

8 
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Type XII 


25 

"2 

(7 - 

■ 7^) ( 7 '^ 

“ 7 ' 

7 /I- 

- 7 

" yr 

“ 7 '^ 

cos (3f 

+ 3.J + 3f' + 3a/) 

• ¥ 

(7 - 

7^y ( 7 '" 

- 7 ' 

®) cos 

(3f 

+ 3<;^ 

+ 3f ' 

+ 3a/ + 

0) 

8 

(7^ “ 

' 7 ®) ( 7 ' ■ 

- 7'^ 

n2 

) cos 

(3f 

. 3a) 

+ 3f' 

+ 3aJ ' “ 

0) 

15 

4 

(y - 

2y3 + yS-j 

y.5 

yr- 7 ^ 

Vi 

- 7 ' 

^ cos 

(3f + 3.a 

+ 3f' + 3fj' + 20) 

15 

4 

7 ^( 7 ' 

- 2y'3 + 

7’^) 

v7- 

y2 

/r- 

7 ' ^ cos ( 3f + 

3 ..) + 3f' + 3a/ - 26) 


+ ( 1 - y' ^ cos (3f + 3(^j + 3f' + + 3f ') 

n 

+ 7^(1 “ y ' cos (3f + 3a; + 3f' + 3a/ - 37) ; 

o 

Type XIII 

+ “ (7^ - 7^) (y'^ - 7 '^) v'^1 ^ 7^ - 3f' - 3a/) 

+ -— ( 7 ^ “ 2y^ + 7 ^) {y‘^ “ 27 "^ + 7*^) cos (3f + 3a; - 3f^ “ 3^o' ' ) 

8 

^ ( “ 7^) (y'^ ~ y ' cos (3f + 3oj - 3f' “ 3a/ - /) 

o 

+ (y - 27 ^ + 7 ^) ( 7 ' - 27 ^ + - 7 ^ /T- 7 ^^ cos (3f + 3a> - 3f^ - 3a)' 

+ “■ 7 ^ 7 " 5 yj I -72 1 - y'2 ^ 3 f 4- 3 ^^^ - 3f' “ 3a>' - 27) 

^ ^ (1 “ 7^)^ (1 “ 7 '^)^ cos f3f + 3a> - 3f' “ 3a/ + 37) 

n 

t 7 y'’y'^ cos f3f + 3 o! - 3f' - 3 ./ - 30) . 

o 

Then, by using Cayley's notations, we obtain 


— ) cos (qf + i) 


OD 


OS + sin”! cos (ig + a) , 
p , q 


Ip' * 

— / cos (q'f' + a') = 2 fcos + sit^l 

V ^ / / I p' , q ' 


cos(ig +a) 


Evidently 


7 )" (/)'' (qf + q'f ' + a) 


+ 03 +C0 

EL 

- Of) * OL 


i i ^ 

[cos + sin] [cos + sin] , , cos ( ig + i'g' 


p , q 


(xin) 



cos (3f - 3f ' + a) . 


a) • (7) 
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The coefficients [cos + sm]^ ^ and [cos + sin] , , are functions of the eccentricities 

and can be taken, in a general case, from Cayley’s'tables (Reference 5). 

For the development of 



cos (qf + q'f' + 


, 


12 


On the top, the factors depending on y and y' are given. Each factor is a 
product of a polynomial in y and a polynomial in These factors are 
analogous to the polynomials of Tisserand from the planetary theories On 
the right are the angles that correspond to the y , y> factors and that 
must be added to the arguments. Each coefficient above must be multiplied 
by each coefficient below, and the arguments are added. The common fac- 
ors, m'a /a' for Types I-V and m'aVa"* for Types VI -XIII, are omitted. 


As an example, from the arrangement for Type 


(below) we take, to order and e'2, 


3 C y - 2 V 


') (y' - 2y'3) /T 


i'3 


4 (~T 2g + 2g' 


''0 


There is no necessity to perform the actual multiplication; 
each case. 


it can be done separately in 


Type aO, 

Factors 

+ I (1 - 6>-2 + 6y'») (1 - 6r'2 + ey'^) 
+ 3 (y2 - ^-^-2 _ ^^-4) 

+ 3 (y - 2y3) (y' - 2 y' 3) 


COS terms 

2 t 0 

2t ^ 20 
2v 0 






+ |-e'2 

- 2g' 



- r' 


(1 - e'2)-3/2 

O.r' 


^ 3 , 27 , , 



+ 16^ + ••• 

+ g’ 


■■■ 

+ 2g' 

- ^ e ^ i e ^ + • • • 



4 12 


- 2g + 2 t 

1 - 




N3| cn 


Type 3^2 


cos terms 
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Type an^ 

Factors cos terms 


+ |- (1 - 6 7^ 

+ ey**) (r'2 - r"*) 

2r' 

= 2^' 


+ y (r^ - r"*) 

r'^ 

2t' 

r: 2o)' 

+ 26> 

+ y - r'*) 

(1 - r'^)^ 

2 t ' 

= 2oi' 

- 20 

+ 3 (r - 2r^) 

\/l - y^ y/l - y*^ 

2t' 

2w' 

+ <9 

- 3 (r - 2y3) 

(y‘ - r'^)\/i - - y'^ 

2t' 

^ 2o)' 

- 




Factors 


cos terms 


Type aQj 


+ -L Cy2 _ ^ 4 ^ (.^-2 _ ^- 4 ) 

+ -f- (1 " r"* 

+ -f- y'’ (1 - y' 2)2 

- 3 (7 - 7 ^) 7 ' ^ 7 l - 7 ^ \/f “ 7 '^ 

- 37 ^( 7 ' - 7 '^) >/T - 7 ^ yi - 7 '^ 


2t ^ 2oj 

2r = 2co + 26* 

2 r ^ 2 ^xJ - 26 ' 

2 r = 2 oj + 6 ^ 

2 t 2 oj - 6 ' 



“ g' + 

0 . g ' + 2f^/j 

+ g' + 2f/J 

+ 2 g' + 2 cu‘ 
+ 3g' +2^4^' 

+ 4g' + 2^4 j' 


~ g + 2r 
O.g + 2t 
+ g + 2r 
+ 2g + 2r 
+ 3g + 2t 
+ 4g + 2 t 



x|u» 


Type aOg 


COS terms 


Factors 

_ 5 y 3 + 5 ^ 5 ) - 5 y '3 + Sy- 5 ) 7 l - y 2 /T - 

+ — (1 -7^) (1 - 107^ + 157"*) (67'^ - 207'^ + i 5 r'®) 

8 

(67^ - 207^ + 157®) (1 - 7'^) (1 - lor'^ + isr"*) 

+ ^ (7 - 47 ^ + 37®) (27'^ - 37'®) \/l - 7 ^ /l”- 
+ ^ (27^ - 37^) (7' - 47'^ + 37'®) \/l - 7 ^ /l - 7'^ 

+ ^ ( 7 ^ - 27 '* + 7®) (7"^ - r'®) 

+ ^ ( 7 '’ - 7®) (r'^ - 27 "* + 7 '®) 


T = <X> 


r = o> + 


r = oj “ 


T = OJ + 26* 


T = OJ - 20 


T = o) + 36^ 


T = rxJ “ 3^? 







“ g' 

+ 


+ e' + • • • 

O.g' 

+ Cj' 


+ 1 + 2e' ^ + ■ - • 

+ g' 

+ oj" 


+ 3 e ' + • * ■ 

+ 2g^ 

+ 



+ 3 g' 

+ CJ'' 









“ g 

+ ■ T 

5 15 3 . ,, 

~ ~2 ^ ~-^e ( exact) 


O.g 

+ T 

+ 1 + 2e 2 + • ■ • 


g 

+ T 

--|e + ••• 


+ 2g 

+ T 



+ 3 g 

+ r 



oo|oj 


Factors 


cos terms 


Type aO^ 

+ (7 - Sy ^ + 57^) (7' - 57' ^ + 5 y ' \/l “ 7^ y 1 - 7^2 

+ J- (1 - 7 ^) (i _ I0y2 + lSy 4 ) _y.2) _ 20^.2 + lSy' 4 ) 7. _ ^ + g 

+ -| ( 6 y^ - 20y* + ISy®) (6y'^‘ - 20y'^ + 15y'®) t = co - 0 

+ -^ (7 - 4 y^ + 37®) (7' - 47'^ + 37'*) v/ 1 - 7^ ■ J \. - y '^ T - CO ^ 29 

+ ^ (27'* - 37®) (27'^ “ 37'®) \Jl - y'^ \/l - 7' 2 T - CO - 26 

+ ^ (7^ “ 27^ + 7®) (7'^ - 27"* + 7'®) T = O) + 30 

+ ^ ( 7 '* - 7®) ( 7 "* - 7 '®) T = CO - 30 



Type afig 


Factors 

r (7 -* + 5y^) ( 7 '^ " 7 '^) \/ 1 " Vl " 7' ^ 

^ (1 - y2) (1 - 10y2 + 15^4) (^.4 _ ^.6) 

^ (67^ - 207'’ + 157 ®) ( y ' - y '^)^ 

- (7 - 47^ + 37®) 7'^ •Jl - y ’^ \/l - 7' ^ 

^ (27^ - 37®) r'(l - 7'^)^ Vl - 7 ^ \/l - 7'^ 

^ (7 - 7 ®)^ 7'® 
r (74 - 7®) (1 “ 7' 


COS terms 


T - CO 


r. - CO 0 


r - CO - 0 


T ^ CO 26 


r ~ CO - 26 


T ~ CO 3*9 


T - CO - Z6 


0 (exact) 



- e ' + • • ‘ 

+ 1 - 6e ' 2 

+ Se' + • • ■ 


+ 



0 . g ^ + 3oj 

• • '*■ g ^ 

+ 2g^ + 3o)‘ 

+ • • * + 3g ^ + 3 <)d 

+ 4g" + 3co' 

• * ' + 5g ' + 3oj 


(exact ) 


- g + r 
O.g + r 

+ g + T 

+ 2g + T 


+ 3g + T 
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Type aQg 


Factors 


cos terms 


+ {y - 57^ + 57^) ( y ' ^ - 7'^) - y '^ >/ 1 “ ? ' ^ 


+ ^ (1 - (1 - lOy^ + 157) (r' " y'^)^ 


r - 0 } ^ 0 


+ ^ (67 - 207 + 15 y®) (y"* - r'®) 


T - (jO - 6 


- ^ (y - 4y3 + 3y^) (y' - 2y'^ + y'®) / 1 - y^ '/l “ y'’ 


T = + 26 ' 


^ (27^ •“ 37^) 7^^ Vl - 7^ \/l - 


r = c<j - 2£? 


^ y2(l - y2)2 (1 - y'2)3 


T = o) + 3<9 


+ ^ (y"* - y®) y' ® 


T = c/j “ 3<9 




+ 5 e' + 


+ 1 - 6 e '2 + 


- e' + 




0 (exact) 


5 15 3 / 

y e - e f exact) 


O.g + T 


1 + 2 e^ + 


+ g + T 


2 g + T 




+ 3 g + T 



Factors 


cos terms 


T - 

T = 3cO + 0 

r - 3co. - 0 

T - 3co + 29 

T - 3 ^^ “ 29 

T = 3w + 30 

T - 3co - 39 





Type aQjj 


Factors 


y5) 


+ 


(y - y3y2 (1 _ y^2^ _ jO^/2 + 15y'^) 

(7^* - y6) (5^1 2 _ 20y'^ + 157^^) 

(7 - 27^ + 7 ^) (7' “ 47'^ + 37'^) y 1 ” 7^ \/ 1 - 7 ' ^ 


fl - y2)3 (^< - ^'3^2 


COS terms 

T - 3 co 

r = 3<^ + ^ 
T = 3o> - 6^ 

T = 3 cj + 26 ? 
T = 3oj - 26? 
T = 3 cd + 3 ^ 


Ar"* - 7 '®) 


T = 3cj - 36? 



‘ 3 g' “ cu' 

‘ 2 g' - 

- g' “ w' 
O.g^ - 

+ - a;' 


0. g + T 

+ g + T 

+ 2g + r 
+ 3g + T 
+ 4g + T 

+ 5g + r 




Factors 


Type ani3 
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However, in the case of the existing artificial satellites only the long-period terms 
are kept. The first power of e' is kept only in the second Legendre polynomial, because 
of the presence of the additional power of a/a' in the parallactic term. The accuracy of 
the coefficients of cos (O.g + 2 t) makes the resulting development of the disturbing func- 
tion valid for all eccentricities and all inclinations, provided the ratio a/a' is small 
enough to secure fast convergence of the development of Equation 2. Thus we have the 
following: 

Type a [n] j 

Factors cos terms 

+ ^(1 - 6y2 + 6 y *) (1 - 6y'2 + 67"*) ^1 ’ 2 t = 0 

+ 3 ( y ^ - y ^) ( y '^ ~ y'^) ’ 2 t = 2^^ 

+ 3 (y - 2y3) (y' - 2y'^) y/ 1 - y^ 2 t = 0 

Multiplied by 

3 

*7 ”2 ( ^ ^ 

+ e ' cos (2 t -g') + (1-e'^) cos 2 t ^ ~2 ^ ^ * 


Type a [n] 2 

Factors 

+ -| (y^ - y“*) (y'^ - y"*) ' “f , 

+ -|- (1 - y2)^ (1 - y'^)^ ■ ’ 

3 4 / 4 5 2 

+ 3 (7 - 7^) (y‘ - y’ y 1 “ 7 ^ \/ 1 “ y' ^ ^ ~ , 

+ 373 3 y 1 - y2 y 1 - y ' 2 . A ^ 

Multiplied by 
7 

+ -^ e ' cos (2t - 3 g' - 2 ^>j' ) + cos ( 2 r “ 2 g' ” 2 ^^) 

-ye' cos (2r - g' - 2c*j' ) . 


COS terms 


2r = 2^/J 


2r = 2o) + 26' 


2r = 2^^ - 20 


2r = 2w + 6^ 


2r = 2ai - 6 


^The second term becomes a secular one for 2r = 0 . 
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Type o [n] 3 

Factors 

+ y (r^ “ y"*) (1 " 6y'^ + 67"*) • y , 

+ y (1 - y2)2 (y'2 _ ^-4^ • y e2 . 

+ y y"* (y' ^ - y"*) ' "f ■ 

- 3 (y - y^) (y' - 2y'^) 7 1 - y^ 7 1 " y'^ • y , 

+ 3 y 3 (y' - 2y'3) ^1 - y2 7l - y'2 • ye^ , 

Multiplied by 

3 , 3 

+ -^ e cos ( 2 t “ g ' ) + cos 2 t + -^ e ' cos ( 2r + g') . 


COS terms 


2r - 2 oo 


2r = 2oj + 2^-’' 


2r = 2^J “ 20 


2r - 2o) 0 


2 t = 2 oJ - 


a[a]^ 

Factors 

+ y (1 - 6y2 + 67^) (r'2 - ^' 4 ^ . ^ f • 

+ -| (^2 - 7^*) 7^^ * ^ Y « 

+ y(y2 - ^ 4 ) +ye2^ , 

+ 3 (7 “ 27 ^) 7 '^ y 1 - 7 ^ yi “ 7 '^ * Y ’ 

- 3 (7 - 27^) (7' “ 7'^) v/l - 7^ * >/l - 7'2 . +Ae-j , 

Multiplied by 


cos terms 


2r - 2<y 


2r = 2co' 4 27^ 


2t 2'^y “ 2f^ 


2t = 2<'*j' + O' 
2t ^ 2<y - 0 


1 7 

-ye' cos (2t + g') + cos (2t + 2g') + y e' cos (2t *• 3g') . 



27 


Type a [n] 5 


Factors 


+ |-(r^ - r'‘) (y''^ - r"*) ’ 

7 . 0 ^ 2,4 52 

+ f (1 - y"') y '* • -y e , 


+ 3 ^4 (1 - ^- 2^2 . _|_e 2 , 




3 (7 


- 7 ^) r ' ^ yi “ • /i ' "X' ^ * T ^ 


- 3y^ (y' - y'^) - y^" • ■ 2 

Multiplied by 


cos terms 

2t = + 2"' 

2t = 2o) + 20 + 2^^' 
2 t = 2to - 20 + 2 ^' 

2t == + ^ + 2<^* 

2t = 2oJ - /9 + 2t*j' 


- ^e' cos (2t + g') + cos (2t + 2g') + ye' cos (2r + 3g') 


Type a[0]g 
Factors 

_ I - 5^3 + 5^5^ - 57' 3 + 57' ') ■ (- - T ’ 

+ 4 - (1 " y^) ("l “ ^^y^ ^ 157 '*) ( 67 '^ - 207 ''* + ISy’®) • i^~~2 ^ ~ ~8 ^ ) ’ 

o 

+ I ( 67 " - 207 “ + 157®) (1 - y'") (1 - loy'^ + 1 S 7 "*) • (- |« - T ' 

+ ^ (7 - 47 ^ + 37 ^) (2y‘^ - 37 ' y 1 - 7^ \/l ~ y'^ ■ ^ ^ ^ ) ■ 

+ (27^ - 3/5) ( /' - 47'5 + 37'5) 7 • " 7^ v/l " 7'"' ' \ ^ ~ ^ ' 

+ ^ (7^ - 2 /'* + 7 ®) (7"* - 7'®) • I" ® ’ 

+ Y - 7 ®) (7'^ - 27"^ +7'®) • (--f e - ^ e 3 ^ , 

Multiplied by 


cos terms 

T - CO 
T - CO + 6 
T - CO — 6 

T = OJ + 20 

r ^ CO - 20 

T ~ CO + 30 

T - CO ^ 30 


+ COS ( T + g ' + oj' ) . 
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Type a [0] ^ 

Factors cos terms 

+ I c> - 57^ + 5^5) _ 5^, 3 + 5^, 5^ 7i - y2 yi - y'2 . 5. e _ ^ ^ ^ ^ 

+ I (1 - r") (1 - 107^ + ISr'’) (1 - 7'^) (1 - 107'^ + 157'“) • (-|e -f e^), r = a, + 0 

+ -f (67^ - 207“* + 157®) (67'^ - 207'“ + 157'®) T ^ <0 - 0 

+ T (7 - 47^ + 37®) (7' - 47'® + 37'®) a/I - \/l - y '2 • (-y e - . r = + 261 

+ ^ ( 27 ® - 37®) (27'® - 37'®) n/i - 7 ® v/l - 7'® • (-y'^ - , r = CO - 20 

(72 _ 2^4 , ^6) ^^,2 _ 2^,4 , ^.6^ . ^ r = .7 30 

+ ^ (y^ ~ y^^ (y’^ - . /_ 5 15 s\ 

8 ^^ ^ ^\T T /’ T ^ G) - 30 

Multiplied by 

+ cos (r - g ' - ). 


Type a[0]3 

Factors 


cos terms 


- (y - Sy^ ^ Sy^) (y'^ - y'^) /l - y2 ^ 2 . A ,4 > 15 ^3 

\ 2 8 

+ ^ (1 - 7 ^) (1 - 107® + 157“) (7'“ - 7'®) • (- I - y ^ . 

+ ^ (67® - 207“ + 157 ®) (7' - 7'®)^ • (- y e - ^ , 

+ X ^ /l - 7 ® /l - 7'®'’ • (- y e - ^ e-j , 

- ^ ( 27 ® - 37®) 7'(1 - 7'®)' /T^ /l^V^ ■ ("T^ " T 

* ^ (y ~ 7®)^ 7 '® • (- I e - y e3^ , 

+ y (7“ - 7 ®) (1 - 7'®)^ ■ (■ I ■ T • 

Multiplied by 


T - CO 


T - CO ■{■ 0 


r - o) - 0 


T = CO + 20 


T ~ 00 - 20 


r = o) + 30 


r - CO - 30 


+ cos (r + 3g' + 3(^' ). 
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Type a [Q] g 

Factors cos terms 

+ ^ (7 - 57 ^ + 57 ®) ( 7 ' ^ - y' -Jt - 7^ v/l " y'^ ■ (“ T ® ’ T - a. 

, ^ (1 - ^2^ (1 - 107^ . 157 ^) (7' - r'^)^ • (- I e - f e 3 ) . r = . . ^ 

+ y ( 67 ^ - 207 “* + 157 ®) ( 7 "* - r'®) • (- I ^ ' T - CO e 

- IL ^ y -^ y ^ , 37®) (7' - 27'^ + r'®) yr^ (-| - - f 

+ ^ (27® “ 37^) 7' ® y 1 ” 7^ y 1 ” 7' ^ y ^ y ’ T ti) 2 

r = „ , 3< 

.f (r^-7)r'« - (-|.-f .’)■ r = --3< 

Multiplied by 

+ cos (r - 3g' •“ ) . 


r>;p^ a[n]jQ 

Factors 

_ i|- (y 3 - 75 ) (y' - 57'® + 57'®) ^f\ - 7 ^ yi - 7 '^ ■ (- ^ e®) • 

+ ^ (7 - 7®)^ (67'^ - 207"' + 157'®) • (- ^ e®)- 

+ X ~ ^ ' (~ ^ ’ 

- ^ (7 - 27® + 7®) (27'® - 37'®) yi - 7^ yi " 7' ^ • (- ’ 

+ ^ 7 ®( 7 ' - 47 '® + 37 '®) yi - 7 ^ \/i - 7 '^ ■ (- y . 

+ ^ (1 - 7^)^ (7"* - 7 '®) • (- y 

7 ®( 7 - - 7 '®)^ • (-f e®), 

Multiplied by 


cos terms 

T - 3 co 

r = 3 <^a) + & 

T - 3 i^ “ 0 
T - ^oj 26 
T - ~ 20 

' 7 ' = 3 oj + 36 ^ 
r = 3<^ ■“ 36 


+ cos (t + g' + O)' ) . 
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Type a[n] 


11 


Factors 


15 


(r^ “ 7^) (y' - Sy‘^ + 5r'^) \/i - y^ /T 


ff-’). 


cos terms 


T ~ 3(0 


15 

8 

(7 - y^O 

' (1 - 7'^) (1 - 

lOy'^ + 157 ' 

*> • (- f 


T 

= 3^ 

+ 6 ;^ 

15 

8 

(74 _ 

( 67 ' 2 - 20y'* 

+ 157 '®) ■ (- 

f =>). 


T 

2 ^ 3(^1 

- ^7 

15 

4 

(7 - 2y^ 

+ yS^ _ 4 ,^, 

^ + 3y'®) /T 

- 7 ^ \/l - y'^ ■ 


r 

— 3'^ 

+ 26< 

15 

4 

7®(2y'3 - 

3y'®) v/l - y2 

/r^'3 . 

(-f'“)' 


r 

= 3'^2J 

- 2 // 

15 

8 

(1 - 

( 7 ' - y'^)^ • 




T 

= 3 2 J 

+ 3^2 

15 

8 

7®(7'^ - 

y'®) • (- f e3 

). 



T 

= 3<'2.‘ 

- 3^:^ 


Multiplied by 

+ cos (r - g' - ) . 


Type a[n]j 2 

Factors 

- y (7^ - 7^) iy'^ - y'^) -JY-T^ yi - r'2 ■ ^ ‘ , 

+ ^ (7 - 7^)^ (y"* - 7'®) ■ (- -^1 , 

+ ^ ( 7 “ - 7 ®) ( 7 ' - 7 ' 3)" ■ (- f , 

- (7 - 2y3 + ^5) ,^-s ^ 

- ^ 7®(7' - 2y-3 + y'S) /l 2 /T-Y^ ■ ^ , 

+ I (1 - y2)3 ^,6 . 35 ^ 

^ I ^-^n- X'V •(- -]? o 3 ), 

Multiplied by 


cos terms 

T :z 3co 
= 3^iJ + 

T - 3oj ~ 

T ^ 3o> + 26^ 

r = 3^0 ~ 20 
T = 3x.; + 3 /:^ 
T - 3^2) “ 30 


+ cos (7 + 3 g' + 3 i)') . 
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Type a [^^1 1 3 
Factors 

+ ^ (^3 _ ^5^ (^-3 _ ^.5^ /T^ • (- f «") - 

, 21 (^2 _ 2^4 , ^6) (^.2 _ 2^4 , ^.6) . (-f e3). 

o 

+ ^ (r-* - y®) (r'" " r'") • (" f^")- 

+ ^ (r - 2y3 + yS) (^'-2y'3 + y' 5) /T^ ' (" 

, ^ ^5^.5 • (- f • 

+ -g- (1 - (1 - y’^)^ ■ (" X ■ 

^|yV^ • (-f e^). 

Multiplied by 


35 

8 



+ cos (r - 3 g' - 3 c*j' ) • 


COS terms 

r - '^Oi^ 

T = 3oJ + (9 

T = 3<>2 “ 0 
T = 3^^ + 2^ 
T - 3<^ “ 2^ 

T = 3oj + 36> 

T = 3^*2 "" 3^ 


The standard equations for variations of elliptic elements may be used in connection with 
this development. In particular, for the variation of the perigee height the equation is 


d Sq x/l - 3 aO 

dt na^e 

The proposed scheme leads to a form of the development which is standard m celes- 
tial mechanics. It differs from Kozai's scheme only by arrangement. The foUowing 
method was used to obtain the development on the automatic computing machine: 


o i 2 

G = cos^- ^^^ 2 

• 9 ^ 2 ^ 

+ Sin-^ -7T COS^ ~7i 


9 1 ' 2 ^ 

+ cos^ sin^ 


. o 1 • 2 A 

+ sin^ ^ Sin^ -;r 


+ sin i ’ sin i * 


sin 1 * sin 1 


cos (cj + 0 “g* ~ O)* ) 

cos (co ~ ^ + g^ + ) 

cos (o) + ^ + g^ + Ca)* ) 

cos (co - 6 - g' - CO*) 

cos (co - g’ - CO* ) 
cos (co ^ g* + CO* ) , 
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Then, 


H - G . 2 ^e' sing' = G . -g-). 


(2 h 2 - + 3e' cos g') - J-h) (1 + 4e' cos g') 


The terms of K are then sorted into four groups according to the coefficient of a, in their 
arguments: the 2 a. terms are multiplied by (+ 1 e^); the 0 • a, terms are multiplied by 


/ 3 \ * \ 2 tcxuio <ixt? mumpiiea Dy 

^ i i y - terms are multiplied by (- |.e - ^5 ^nd the 3 a. terms are 

multiplied «/ 

In order to include the main effect of the lunar and solar perturbations in the method 
used at the NASA Space Computing Center (Reference 6), the development of the disturb- 
ing function must be obtained in terms of the eccentric anomaly ^ of the satellite. By 
eeping the second harmonic only and using the formula 

(f) cos (2f + 2r) = (1 + ,32)-' cos (- 2i + 2r) 

“ 4/3^ cos + 2r) + 6/5^ cos 2r 


4/3 cos (+ c + 2r) + cos (+ 2^ + 2r) 


the following computational scheme is obtained, 
ties and all inclinations. The common factor (m^ 


This scheme is valid for all eccentric i- 

a^a'^)(l 4 /3^) ^ is omitted. 


Type 

Factors 

+ ^ (1 - 6y2 + ~ ^ 


cos terms 
2t = 0 






Nd|OJ K)|U) ro|o.i 
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Type aflj 


Factors 

C05 terms 

( 7 ^ - 7 *) (1 “ 67 '^ + 67 '*) , 

2r = 2 o) 

(1 - 7 ^)^ (r'^ - 7 "*) . 

2r = 2c»J + 2^ 

y4 ^y'2 _ ^-4) _ 

2r = 2oj - 2^ 

(7 - 7^) ( 7 ' - 27 '^) yi - 7^ \/i - r'^ . 

2t = 2co + (9 

(r' “ 27 '^) x/l - yjl - 7^2 ^ 

2r = 2o> - ^ 



. 3 r 



4-e' . ••■ 

‘ g' 


1 + 1 e'2 + ... 

O.g' 


3 , 



+ _e' + ... 

+ g' 

+ 


- 2f + 2 t 




“ 4/3 


“ ^ + 2t 

2 



+ 6/3 , 


0.^ + 2t 

- 4/3 


+ ^ + 2t 

+ 1 


+ 2^ + 2t 



Type aO^ 


Factors 

COS terms 

1 (1 - 67 ^ + 67 '*) ( 7 '^ - r"*) . 

2t = 2o' 

I ( 7 " - r") r"* , 


2t - 2oo’ +4 

1 (^2 _ ^4) 


2r = 2<^' “ - 

3 (7 “ 27 ^) y'^ Jl 

- 7^ Vl - 7 ' ^ . 

2t ^ + t 

3 (7 _ 2 y 3 ) ( 7 ' - 7 

v/l - 7^ x/l - 7 '^ . 

2t = 2(^' - < 


0 f exac f ) 

O.g' 


1 ' + 

“ — e + * * * 

+ g' 


2 


1 _ i e<2 + ... 

1 2 

+ 2g' 


t y e' 4 ■■• 

+ 3g' 




2 

+ /i 


- 2^ + 2t 

1 

K5 

1 

K5 


- ^ + 2t 

2 4 


O.f + 2r 

1 + 4/5 + p 


- 2fi - 2/3^ 


+ ^ + 2r 

+ 


+ 2^ + 2r 



Factors 


cos terms 


Type aOg 


- y4-) 

2 _ , 

2t - 2o) 2^' 

y2)2 ,^,4 

- 

2r =2 2^^ + 20 -h 2<^‘ 

- y'2)^ 

- 

2r = 2oj - 20 ^ 2rj' 

"X^ ^ 

/T - 7^ \/l - 7'^ , 

2r = 2oj 0 + 2co' 

- r'^) 

\fl - yi - y'2 ^ 

2r ' 2oj “ 0 + 2(^' 


0 f ejcac f ^ 

O.g' 


1 , 



- 2 e' + •-. 

+ g' 


1 - 1 e'2 4 ... 

+ 2g' 


. 7 , 



2 '^ *** 

+ 3g' 


- 2^ + 2r 
“ ^ + 2t 
0.^ + 2t 
+ ^ + 2r 


+ 2^ 4 2t 
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THE INTEGRATION PROBLEM 


The method used at the NASA Space Computing Center is based on Hansen’s idea of 
separating the perturbations in the orbit plane from the perturbations of the orbit plane. 
A typical differential equation of the method takes the form 

dQ cos 

TF ^ ^ ' 

ci E / , s 1 n 

a - iE + jtd+kr2+ + jVx/ + + pF , 

where f is considered as a constant and 

(V = (E - Eq) , jx' = + n' (t - t^) , 

Q ^ (E - E„) , a;' r -i (t - t^) , 

i)' . -U; 4 nVO; (t - t,) . 

The time t can be eliminated by means of Kepler’s equation, 

E-esinE 

and, after several transformations involving Bessel functions, the right-hand side can be 
reduced to an integrable form. However, instead of this complicated method, the integra- 
tion can be performed in such a way that E and g' are both kept in the argument. A 
similar method was used by Hansen (Reference 7) in his development of the perturbations 
of Encke’s comet by Saturn. 

From Kepler’s equation and from 


r' Ro + ~ ^o) 


it can be deduced that 


d f 


m ( 1 - c cos i;;) 


^ who re m 



and 


dQ c)Q r)Q dg' f)Q 3Q 

77” ~77^ ^ 7T”^r * ~T” '7~ ^ ^ ' — 7 (1 - e cos 

d(^" c)7- r)g cit ()<*; ()g ^ 


Thus the problem is reduced to the integration of the linear partial differential equation 


c^Q c)Q 

— 4 m ( 1 - e COS ) 

Og' 



cos 

A 

s i n 


a 


X . 


— — - ™ ^ tMAVA f T iJULW Ai 

extended their calculations to the case of the paddle-wheel satellite, Explorer VI 
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Figure 1 - Solar and lunar perturbations on the perigee height of 
Vanguard I (1958 /i^) 

(Reference 1). Explorer VI has an apogee of 48,700 kilometers, a perigee of 6640 kilo- 
meters, and an orbital inclination of 47.3 degrees; it was found that this highly eccentric 
orbit produces substantial lunar and solar perturbations which decrease the perigee 
altitude rapidly and shorten the orbital lifetime from several decades to a probable value 
of two years. 

The interesting work of Kozai and Whitney encouraged further exploration of the pos- 
sible lunar and solar effects on perigee height for satellite orbits of large eccentricity. 

In general, both the eccentricity and the perigee height vary with time as a result of these 
effects. The amplitudes, frequencies, and relative phases of the variations are deter- 
mined by the orbit parameters, among which the hour of launch is of considerable 
importance. For a special set of launch conditions, and for representative orbit parame- 
ters, the perigee height may he made to rise steadily over the course of several years at 
a rate of 1 kilometer per day. Thus, the sun and the moon may provide a substantial 
perigee boost for the satellite under properly chosen circumstances. For other condi- 
tions the perturbations may be minimized to obtain a relatively stable orbit. These con- 
siderations may be of importance in deciding the launch programs for future satellites 
with highly eccentric orbits. 



PERIGEE HEIGHT (KM) 


6Y 
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From the perturbing function developed herein, the rate of change of perigee height 
was found by the method of variation of constants. Letting q be the perigee height, 

dq (1 - e)^ ^0 dQ 

dt nae 3 ct nae 

where cj is essentially the mean anomaly g of the satellite, rx; is the argument of 
perigee of the satellite, and n is the perturbing function. 

The first-order perturbations of the gravitational effect of the sun and moon were 
developed, and several assumptions which considerably simplified the expression for the 
potential were made. It was assumed that the moon’s orbital plane is coincident with the 
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Figure 2 - Solar and lunar perturbations on the perigee height of 
Explorer VI (1959 
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sun’s. Only the long- duration effects on perigee altitude were considered; this eliminated 
terms with periods of the mean anomaly (and fractional parts thereof) of the satellite and, 
for the most part, of the moon (i.e., g and g^). 

Hence, the expression dq/dt becomes 

dq _ ^ 

dt nae ' 

which is a trigonometric series whose arguments depend on various combinations of 2 oj, 

^0 > ^ terms of these elements several possible resonance effects can be 

recognized, two of which are associated with the following conditions: 


2^-2 

- 0) = 0 , 

( 10 a) 

2co + 2 ( A-g 

- 0) = 0 . 

( 10 b) 


where lo ^ \ , e are the average angular velocities of a;, , e , respectively. 

These resonance conditions have a simple interpretation. For example, in Equation 
10 a, ^0 - 0 represents the longitude of the sun relative 1 o the line of nodes; and a. , the 
position of the perigee in the orbital plane, is also defined relative to the line of nodes. 
Therefore, in a system in which the line of nodes is fixed, the satisfaction of the reso- 
nance condition (case 10 a) signifies that the mean angular velocities of the sun and 
perigee are equal; that is, the line of apsides follows the sun. In this case the orbital 
perturbations produced by the sun are clearly maximized. In Equation 10b, the sun and 
the line of apsides have the same period of revolution but in opposite directions. Again, 
it is clear that the solar perturbations will be maximized. 

It was found that the effects on perigee height may be maximized or minimized by 
choosing suitable values of the orbital inclination and the time of launch. A long-period 
effect occurs when the inclination to the equator is near 63.4 degrees — the critical angle 
at which the motion of the argument of perigee is small. There is a 2 co term in dq/dt 
which at this inclination produces a near-resonance effect, causing the perigee height to 
change almost secularly. For an orbit with an apogee of 46,550 kilometers and a perigee 
of 6650 kilometers, the rate of change of perigee is about 1 kilometer per day, as shown 
in Figure 3. The precise magnitude of the rate of change depends on the initial argument 
of perigee. The hour of launch does not affect this result. 

At angles of inclination other than 63.4 degrees, a variety of effects may be obtained 
by a suitable choice of the hour of launch. Selecting the hour of launch is equivalent to 
selecting the longitude of the ascending node, with any value available once in 24 
sidereal hours. The results of three different launch times are shown in Figure 4, 
employing the same apogee and perigee as in Figure 3 with an argument of perigee equal 


G-4 back 
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to 135 degrees and an equatorial inclination of 28 degrees for February 1, 1960. Curve A 
corresponds to a launch time of 7 hours U.T. on February 1, curve B to 23 hours U.T., 
and curve C to 13 hours U.T. Cases A and C show rapid initial variations of perigee 
height — advantageous and disadvantageous. Case B represents a relatively stable orbit. 

Curve C represents the addition of drag to case C. It rises initially above the solar 
and lunar perturbation curve because the drag decreases the period and the eccentricity, 
and these changes in turn decrease the solar and lunar perturbations. It is interesting to 
note that the lifetime for a satellite with the parameters of Figure 4 is 25 years in the 
absence of lunar and solar perturbations and approximately 1 year with these perturbations. 
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